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ABSTRACT 


In  this  paper  tha  non-ayimeetric,  free,  elaatic  vibrations  of  thin  domes 
of  revolution  are  studied.  It  is  aasueed  that  the  frequency  is  low»  The 
asymptotic  approximations  previously  gives  by  the  writer  are  used  to  estimate 
the  general  solution  to  the  shell  vibration  equations  at  low  frequencies. 
Approximations  for  the  low  natural  frequencies  and  nodes  are  derived  systemat¬ 
ically  under  a  variety  of  edge  conditions.  Low  natural  frequencies  are  found 
only  when  the  edge  conditions  impose  no  forces  tangent  to  the  shell  surface, 
-fnen  tha  edge  is  free  (and  only  then)  Rayleigh's  inextensional  frequencies 
are  recovered.  Tor  certain  ether  edge  conditions  rew  natural  frequencies  are 
found  that  are  above  Rayleigh's  frequencies  but  still  low  compared  ».r.  with 
the  lowest  mewbrane  frequency.  The  displacement  modes  associated  with  these 
new  frequencies  are  mostly  of  inextensional  type.  The  general  results  are 
apoliad  to  estimate  these  now  frequencies  for  spherical  doees. 
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1.  INTRODUCTION 


lhe  inextensional  vibrations  of  thin  shells  wert*  first  studied 
by  Rayleigh  iij1,  and  since  that  time  his  procedure  has  often  been 
used  to  estimate  natural  frequencies  for  various  shell  shapes.  The 
frequencies  ootained  by  this  procedure  ara  such  lower  (for  a  thin 
shell)  than  those  predicted  by  any  ether  method  and  are  therefore  of 
great  practical  interest  in  such  applications  as  tents,  parachutes, 
and  metal  or  plastic  containers. 

Kowevsr,  there  are  good  reasons  for  skepticism  conceding  the 
generality  of  Rayleign's  procedure.  Tor  cxnsple.  Love  [2]  has  shown 
trm  tne  modes  satisfy  neither  the  sstion  equations  nor  (with  a  few 
exceptions)  the  edge  conditions.  Also,  Arnold  and  rfarfccrton  [3? 
observed  that  Rayleigh’s  procedure  gave  good  agreement  with  experiments 
in  seme  cas-s,  but  that  changes  in  the  edge  conditions  -ould  cause  enor¬ 
mous  changes  in  the  lowest  measured  frequencies  and  coepletely  destroy 
the  agreement.  It  appears,  therefore,  tnat  we  do  not  under* tend  these 
vibrations  as  well  as  we  ought  to. 

In  the  present  paper  we  shall  show  how  i nex tans ion a 1  modes  may  be 
systematically  derive.)  from  a  general  theory  of  shell  vibration,  without 
explicitly  assuming  that  the  mode  is  inextensionai.  Rather,  it  is  merely 
assumed  that  the  frequency  is  low  (  in  a  sense  that  will  later  be  made 
more  precise),  and  froo  this  assumption  inextensionai  nodes  and  fre¬ 
quencies  are  derived.  This  change  of  procedure  is  Important  for  two 

•Numbers  in  brackets  deaignate  References  at  end  of  paper. 
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re^soni.  rirst,  we  find  that  inextensionau  aod*a  can  be  derived  only 
for  certain  edge  conditions,  and  this  sheds  light  on  the  questions 
raised  in  [2j  and  [3].  Second,  the  new  procedure  leaves  the  w;iy  open 
to  find  all  low  frequencies,  whereas  Rayleigh's  procedure  is  limited 
vj  frequencies  for  which  the  modal  bending  energy  greatly  exceeds  the 
modal  stretching  energy.  For  certain  edge  conditions,  we  shall  find 
inextensionai  modes  with  frequencies  different  from  those  obtained  hy 
Tayloigh. 

To  demonstrate  the  procedure  in  a  context  general  enough  to  he 
convincing  but  simple  enough  to  avoid  unessential  manipulations ,  we 
consider  a  general  dome  of  revolution  executing  snail,  non-syrametric 
vibrations.  We  shall  use  the  approximations  obtained  by  the  author  f«l 
vo  write  down  an  approximate  general  solution  of  the  differential  equa¬ 
tion  system  when  the  frequency  is  low.  This  solution  is  substituted 
into  the  boundary  conditions,  the  resulting  frequency  determinant  is 
solved  and  the  ratios  of  the  arbitrary  constants  are  found.  This  en¬ 
tire  process  is  carried  through  for  four  different  edge  conditions, 
starting  w.th  a  free  edge  and  proceeding  at  each  stage  to  the  "freest" 
of  the  remaining  ^dge  conditions.  The  frequency  increases  with  each 
new  edge  condition  until  we  exhaust  all  edge  conditions  for  which  low 
frequencies  can  be  found. 

for  the  two  "freest"  edge  conditions  this  procedure  gives  com¬ 
plete  estimates  of  the  mode  but  only  an  order-of-aagnitude  estimate  of 
the  frequency.  To  find  frequency  estimates  we  use  Rayleigh  s  Princiole 
for  these  cases.  In  the  other  two  ceses  explicit  estimates  are  found 
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for  the  inextensionai  frequencies.  The  general  formulas  are  applied 
to  a  spherical  dome,  and  numerical  results  are  obtained  for  the  pre- 
viously  unknown  inextensionsl  frequencies. 


2.  FUNDAMENTAL  EQUATIONS  AND  SOLUTIONS 

We  shall  adopt  as  our  starting  point  the  equations  of  thin-shell 
theory  propounded  by  Sanders  [5]  and  modified  by  the  inclusion  of 
translational  (but  not  rotational)  inertia.  All  effects  of  transverse 
shear  and  thickness  change  are  omitted  in  this  theory  w*  may  writ* 
the  system  in  dimensionless  form  as  in  DO. 
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The  motion  equation*  are 

(1  (r«  '  f(n  -  nfl«)  t  Hn  J  * 

SS  SS  0*1  35 

♦  ♦  l/2MBS9(rs“‘  -  p#"» )}  =  0  (7) 

(1  -  v2  )'*{n  *  ♦  2fn  -  tta  }  *  Q2v 

^6  S®  00 

♦  eJ{q9r6'1  ♦  l/2t(r9'»  -  p#‘*  =  C  (8) 

(1  -  v2ri{Qssr,*’1  ♦  n  r  "*)  -  a*w  -  c2{q  '  *  fc  ♦  Mq  )  =  0  (9) 

88  8  *  “  9 

H«r«  in  dimensionless  fore  u,  v,  and  w  are  tit*  meridional, circumferential, 
and  normal  displacements,  y««»  Ya«  and  y  _  are  the  middle  -surface  strains, 

89  39 

b  and  b  the  rotations  and  k  ,  k  and  k  the  curvature  charges.  The 

•*  o  5b  SO 

n*s,  m's,  and  q’s  are  the  direct  (membrane)  stresses,  bending  romeni.s_,  and 
shears,  respectively.  Also  r$  and  rg  are  the  principal  radii  of  curvature, 
v  is  Poisson's  ratio  and 
a  =  uR(p/E)«^ 

t '*  ~  h2/[12R*(l  -  „-£)]  «  1  (13) 

f(o)  *  rfl  1  cot  4,  H(o)  =  ar0"‘  esc  4 
where  u  is  the  frequency,  p  the  mass  density,  E  Young's  modulus,  h  shell 
thickness  (assumed  constant),  R  a  length  characteristic  of  the  radii  of 
curvature,  m  the  circumferential  wave  number  and  4  the  angle  between  che 
normal  to  the  shell  and  the  axial  direction.  Primes  denote  differentiation 
with  respect  to  o  ,  which  is  dimensionless  arc  length  along  a  meridian 
The  boundary  conditions  at  an  edge  have  been  given  by  Sanders  f5) 
and  consist  of  prescribing 
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The  principle  of  conservation  of  energy  for  the  vibrating  shell 
states  that 

%  *  h  -  h  *  0 

where 

Ek  *  fl2/(u2  ♦  v*  ♦  w2)  r^sin  *do  *  B2K 

Es  =  (1  -  va>*V< nssyss  ♦  n#etee  ♦  2nseYs#)  r9  sin*  da 

EB  =  c*^"sskss  *  a99,c9#  2H,s#ks8)  **e  sin*  da 

and  the  integrals  are  extended  over  a  neridian. 

He  shall  now  describe  the  approximations  that  we  shall  use  for 
the  solutions  of  this  system.  The  system  is  linear,  of  eighth  order 
and  has  singularities  where  sin  <  =  0.2  We  limit  ourselves  to  tne 
case  where 

■*€  «  1. 

Tour  of  she  eight  solutions  vary  rapidly  with  o  (i.e.  along  a  neridian) 
and  are  called  bending  solutions,  and  four  vary  srucb  more  slowly  and 
are  called  membrane  solutions.  The  approximations  for  tha  bending 
solutions  are  quite  different  from  the  approximations  for  the  membrane 
solutions.  Two  solutions  of  each  type  are  singular  where  sin  *  *  0. 

We  shall  now  list  the  approximations  to  the  four  bending  solu¬ 
tions,  first  near  sin  *  *  0,  then  for  sin  *  f  0.  The  latter  are 


2  W«  assume  that  sin  *  *  0  at,  end  only  at,  the  axis,  sad  that  the 
apex  of  the  dome  is  of  second  degree. 

S 


(12) 


(13) 


linear  combinations  of  the  approximations  obtained  in  [4]  for  the  case 


3r  <  1. 


For  sin  4  *  0: 
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For  sin#  t  0: 
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(18) 


H  i  [(1  -  n/)cj1/"  (2xr.sin#rl/2(r  ~Z  -  3Z)*3/" 
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G  =  (2  ♦  v )r  ~s  -  r~l 
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In  deriving  these  formulas  we  have  assisted  that  a  is  measured  from 
the  apex  of  the  dome.  Also,  R  is  chosen  as  the  common  value  of  the 
principal  radii  of  curvature  at  the  apex.  Then  rs(0)»  rft(0)  *  l  and 
the  definition  of  x  for  sin  4  *0  is  a  continuation  of  that  for  sin  440. 

These  approximate  forauias  were  obtained  by  an  asymptotic  analysis 
of  the  fundamental  system  of  equations.  The  analysis  indicates,  and 
we  shall  assume  henceforth,  that  the  errors  in  these  approx imat ions 
are  all  0( a1 ).  For  example,  we  could  write 

w  s  AjHe’cia)  [1  ♦  A~1tjw^3  *  A2HeTs(a)  [1  ♦l*1nw*2*3  ♦  ... 

%(j)  *  0(1);  j  «  1,  2,  2,... 

and  similarly  for  all  the  other  variables.  General  expressions  for  the 
corrections  are  not  known ,  nor  do  wa  know  even  the  leading  terms  in 
their  expansions  in  A*1  although  theae  could  prasumaL>.>?  be  found. 

Fcr  the  sake  of  brevity  we  shall  refrain  from  indicating  these  correction 
terns  in  our  equations,  upon  the  understanding  that  they  are  generally 
O(a'J).  However,  thsre  are  several,  rather  important ,  equations  in  which 
the  leading  terns  cancel,  and  in  those  wa  shall  write  down  symbols  ror 
the  correction  terns  even  though  we  don't  know  them. 

In  general  we  assisse  that  tares  which  are  OCA*1)  are  negligible  com¬ 
pared  with  those  that  are  0(1).  This  gives  us  a  rationale  for  deciding 
when  theae  ber.ding  solutions  are  approximately  static,  for,  expanding 


dll  toe  approximations  in  powers  of  Q2,  we  sea  that  the  static  asymptotic 

approximations  are  obtained  when 

i£<  0(a*2>.  (1« 

The  approximations  we  shall  use  for  the  membrane  solutions  are  those 

described  in  [*»].  At  this  point  we  assume  that 

a2  <  o  a-1).  (20) 

Then  the  membrane  solutions  separate  into  two  groups.  Two  (labelled 
5  and  »)  are  approximately  the  static  membrans  solutions.  Two  (labelled 
7  and  8)  are  (for  m  >  2)  approximately  the  inextenslonal  solutions. 

When  (20)  is  satisfied,  all  the  quantities  associated  with  these  four 
solutions  are  approximately  static  (i.e.  independent  off!)  and  0(1) 
except  the  direct  stresses  of  the  inextensional  solutions,  which  are 
ol  tne  form 

n  *  A7[fl2n7(0)  ♦  c2n7(e,J  ♦  A^n^0*  «•  (21) 

wnere  n-j  n7  are  0(1),  independent  of  Q  ,  and  are  found 

oy  solving  the  notion  equations  as  a  non- homogeneous  set  of  three 
equations  in  wnich  the  quantities  multiplying  Q2  and  e2  are  known  in 
advance  from  the  conditions  of  zero  strain. 

Of  the  eight  solutions,  four  are  singular  at  +-  0  and  must  be 
tiiscard-sd  for  a  dome.  We  may  take  two  of  these  as  the  membrane  solution 
numbered  6  and  the  inextensional  solution  numbered  8,  and  we  see  from 
(1**)  that  the  remaining  two  are  the  bending  solutions  numbered  3  and  **. 

Thus  we  oust  take 

Aj  =  A4  =  A6  =  Aq  =  0. 

Tne  approximate  general  solution  can  now  be  written 
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*  *  H*Y{AjC(a)  ♦  A2»(a)}  ♦  Aj*^  ♦  A?»/7* 

u  *  A^G^He*  (-A^sCa*)  ♦  A^cCa*)}  ♦  A^5)  ♦  A?u^7^ 

v  *  A"2MGyHev  (A^U)  -  K^U))  ♦  Asv{5>  ♦  A ?v(7) 

bs  *  AHeY  {-A^Ca*)  -  A2«(a*)}  ♦  A^b^5*  ♦  A7bs*7^ 

b8  =  MHeY  (A^U)  ♦  A2*(a)>  ♦  A<b8t5)  ♦  A?b9(7) 

«  *  A2HeY  {A.o(a)  -  A„c(a)  }  ♦  A,*  ♦  Aj 

ss  1  2  5  sa  ~sm  (22) 

m0e  *  a2vH#Y  {Vu)  •  a2cU)}  *  Asmee(i>>  *  A7"eA<7) 

m  «  AM(i-v)  HeY  (A.c(a*)  ♦  A_s(a*>  )  ♦  A,*  ^  ♦  Aj  ^7* 

80  1  2  3  so  '  S) 

qa  *  A3HeY  (Ais(a*)  -  A2c(a*))  ♦  Asq#^*  ♦  A?qs*7* 

qe  *  A2HHeY  {-AjS(a)  ♦  A2c(a))  ♦  A^q^*  ♦  A?q8*7* 

nss  «  A“ifGnHeY  (AAa(a*)  -  A2c(a*)»  +  Asnss(5)  ♦  J*1*^1*  ’) 

na8  *  G*»H«Y  tAjCia)  ♦  A2»(a)}  ♦  A5ne8^5}  ♦  A7(a2n88*Q)  t  ^ng®**’* 
nS0  -  A^MG^HaT  (A^a*)  -  A2c(a*)>  ♦  A5ns8{  5)  ♦  A7{02ns8(fl)>  i”1*ns8<t  *> 

These  approximations  ara  accurate  when  02  <  0  (X-i ),  a  _>  2  and  sin  4  tQ. 

2  -2 

when  Q  ^  a(A  ),  ail  tbe  quantities  ara  approximately  static  (independent 

of  a)  except  tbe  direct  stresses  of  tbe  inextensional  so'ution.  These 

formulas  for*  tbe  basis  for  our  analysis  of  low  frequencies  and  modes. 
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3.  CALCULATION  OF  LOW  FREQUDICIES  AND  MODES 

In  this  Ssction  ws  shall  put  the  gsnarsl  solution  (22)  into  various 
ssts  of  edge  conditions  and  calculate  the  natural  nodes  and  frequencies. 
The  derivation  will  toe  carried  out  in  detail  for  two  cases  but  results 
will  toe  given  for  tbe  rest.  We  shall  begin  with  the  case  of  a  completely 
free  edge  and  then  consider  succesively  "tighter"  sets  of  edge  conditions 
until  the  frequency  is  increased  above  the  range,  Q2  <_0(V^)  in  which 
(22)  is  applicable. 

We  assume  the  edge  is  at  a  *  and  sin  #(oQ)  WO.  A  new  set  of 
constants,  B^,  j*l,  2,  5,  7  is  introduced,  defined  by 
B1  *  Alh  <o0)«*(ao>.  B2  *  A2H(o0)eY(aO) 


B5  =  A»> 
and  we  also  set 


®7  *  A7  • 


(c)  ,  -2  (  t)  ,  -to  (l) 

X  V  .  A  r9  “ss  s  A  nss 
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A  r 
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In  deducing  tbe  natural  frequency  and  evaluating  the  constants  it  is  to  toe 

understood  that  all  quantities  are  evaluated  at  o  *  oQ. 

Case  (I):  Free  Edge.  nm  *  N  .  =  Q  »  *  0  at  o  »  cr 

9V  S  51  0 

The  four  conditions  are  (keeping  tbe  leading  terms  only) 

n  *  B.A'XjIi*)  -  B„A"if8  c (a*)  ♦  B,n  <5) 

**  1  n  2  n  5  ss 

♦  B7(Q2n>8(a>  ♦  A"**n__tA>)  *  0 


=  A  n 


t*> 
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ss 
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In  In  S  si 

♦  B7(a2n88<a)  ♦  A'4[n  (U 


se 


♦  8»  i  2  0 


se 


Q  =  B  A3s(a*)  -  B  A3c(a*>  ♦  B  [q  l5)  ♦  Ate 
*  i  2  se 


t  B?(qs(7)  ♦  Me_„(7,>  =  0 


se 


(23) 

(2“  ) 

(25) 
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a  *  B  A2s(a)  -  B  A2c(s)  ♦  B  a  B  •  <7)  *  0 

••  1  2  5  *•  7  *• 


g  *  l/2re'2(l  -  v2)(3r9~l  -  r,"1) 

If  we  eliminate  B  using  the  condition  (23),  we  obtain  a  three-aquar-e 

syststB  tout  can  ha  written  in  matrix  fora, 

i\  v  in 2y  ^  t  "1  P"  a  I 


l 

X 

(i^X 

11 

12 

13 

1 

X 

(a^x 

21 

22 

23 

i-l/2 

X  A*1 

32 

(ft2  A”1 

TP 

♦  A'^X  )  B  A-2 

13  1*.  2 

♦  A^X  )  B  « 

23  2<*  S 

“*X  ♦  A'^X  )  B 


_J  u 


X  *  n  *5)  -n  (5}(M/f  i,  X  *-nM'!>V(iG  >,  X  *  -X  s(a) 

12  S»  SS  22  *“  n  32  22 

X  =  a  n(a>  -  n  <a>(M/f  >,  X  *  -n„(a)/(fGn),  X  «  -X  a(a) 

13  S8  aa  23  »•  ”  *  33  23 

V  ■  »„‘u  ♦ 

*„  -  -i  »„“’/«=)  • 

X  1  -■  *7^s(a*), 

3*»  as 

Wnen  B  ia  eliminated,  tba  leading  terms  in  tha  coefficients  of  B  cancel  in 
l  2 

equation!  (2«*)  and  (2S).  The  doe inant  terms  in  these  coefficients  then 

arise  fro*  later  teras  in  the  asymptote  expansions  of  n>s,  ns9  and  qs  for 

the  bending  solutions.  These  are  not  know?.  explicitly,  but  we  know  their 

orders  of  magnitude  and  designate  tbe  unknown  functions  n  and  n  ,  both 

11  21 

of  which  are  0(1). 

The  frequency  is  found  by  annulling  tha  determinant  of  this  system, 
witn  the  result 

xx  -  x  x  t  2^2(n  x  -  n  x  >x 

i<*  22  12  2<*  21  12  11  22  3* 


cPa*  = 


XX  -  X  X 
12  23  13  22 


11 


Tne  ratios  of  tba  coefficients  ars  found  to  be 

a  /»  «  A*22l/2m  *7jc(a*> 

l  7  ** 

8  /*  *  A*221''?«..(7)8(a*>  (28) 

B  /B  »  A*1* 

5  7  5 

where  fi  *  0(1)  is  a  constant  that  may  be  datarainad  from  the  system. 

Tba  denominator  in  tba  frequency  condition  (27)  is 

XX  -  X  X  *  n  <5)n  (a>  -  n  (5>  n  <Q>' 

12  23  13  22  **  »«  se  SS 

and  cannot  vanish  bacauaa  tba  diract  strassaa  aaaociatad  with  the  membrane 
and  inextansional  solutions  aust  ba  linaarly  independent*  Thus  in  the  range 
Q2  <  Oil*1)  tne  frequency  condition  can  ba  satisfied  only  whan 
fl  *  Oj/T2 

Since  A  «  1  r  _1^2  whenever  a2  <  G(l_l),  wa  find  that  vhere  is  only  one 
natural  frequency  for  each  *  >_  2  in  the  range  0 2  <_  Oic1^2 ),  and  it  is  given 
by 

a  3  a  jCFg(o j ).  (*9) 

We  cannot  determine  aT  and  8  because  wa  do  not  know  n  ar.d  n  ,  which 
1  5  ll  21 

arc  found  from  toe  second  terms  in  the  asymptotic  expansions  of  n  „ .  n  and 

’  r  ss*  se 

qs.  hence  (29)  is  not  of  much  practical  value  in  calculating  the  frequency. 

However,  a  first  approxistation  to  the  aoda  is  completely  determined  by 
the  coefficients  obtained  in  (28),  even  though  wa  dc  not  know  §  precisely. 
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W  *  = 

>g  *e/ 

q,  *  A(o)  x*ia  t 

%  •  -X*<0)  *ifl  (i  -  (w/H)}  ♦  q  (’>(o) 

0 


“ss  .  f(o) 

*  A~*(o )x  G_(o) 


*  A  2(c)  xGqCo)  ct.s  { (  -  (*/4 )} 


"ee  " 

r  *  2*/^ 


/Mo)  l 
l  a(J^T / 


•-^riSHfer  “■’•<i,<v 

4  =  2~1/2<x  -  *0J  «  -2-l^X  /®»  (re{0')}  *1/2  ^ 

U  i«  noteworrny  that  all  the  quantities  occurring  in  these  fondle,  for  tne 

MOd*  ar#  StatiC  and  «*y  to  evaluate.  The  displacements  and 

rotations  are  dominated  *-y  the  ioexten.ion.1  solution,  the  membrane  solution 

iS  “Dtir*ly  aad  **  (edge-effect)  solution,  have  a  strong 

influence  on  the  stre.s-like  quantities,  eeking  possible  tb.  satisfaction 
of  ail  boundary  conditions. 

p~c*Sur«  •“*  »*•“•*  °°iy  M  ora.r-.f-^p.i,^.  „tfa4t. 

for  rh.  fr.WTOy,  Oa,  a.Uv«r«)  .«  Mtta.,.  of  ,l»  .oj.  tb.,  i,  b„th 
general  and  tore  collate  than  any  previously  known. 


Case  (II);  a 


*  *  Q  =  2  »  o  at  o  =  o. 


Among  tne  possible  edge  conditions  this  is  the  "freest"  except  for  the 
.we  edge  of  Case  (I).  The  analyst.  strongly  resembles  that  of  Casa  (I), 
and  w.  shall  «*rely  record  the  results.  Only  on.  frequency  is  found  in 


a ii  =  *»(i)  cannot  be  found  explicitly  because  of  cancellation  of  the  leading 


terms,  as  was  true  of  a^.  A  couplers  first-approximation  f'r  the  node  is 
found r 

w  *  w (7Ho), 


u  *  u^7^(o) 


v  *  v^7^(a) 


n6  «  »el7)(o) 

bs  3  xcos  C  ♦  b^Ho) 


EX 

ss 

1 

3  -A(a>x 

m  „ 

v 

ee 

sin  {;  -  («/<*)) 


ms0  »  -(1  -  v)H  la)  xcos  c  ♦  »Se^7)(o) 


q_  3  -A2(o>x  sin  C 


q,  3  A  (a)M(a)x  sin  {<;  -  (*/“)} 

V 


n  „ 

f 

ss 

3  -A”2  lo  )Gn(a  >x 

n 

H 

9- 

nee  * 

-A-1 (o )Gn<o )x  cos  (t  - 

sm  c 


> 


(32) 


wnere  ;  is  defined  as  in  (31),  and 


X  3  H(o)  Alo)  .  ra,/l,(o0^  sinl^2#(o0 ) 

— - D(7,(oJ.^-5 -  - 7 - ~3(7>(oa)e 

H(oq)  Avoq)  r  i<l,(o)  sin1  2f(o) 

The  frequency  is  somewhat  aigher  than  in  Case  (I).  The  modal  displacements 
are  wholly  inextensional,  and  the  stress-like  quantities  are  almost  entirely 
derived  from  the  Bending  solutions. 

We  see  that  Cases  (I)  and  (II)  are  quite  similar.  In  neither  case  car. 
we  calculoc*.  the  frequency  directly.  But  in  both  cases  we  have  very  good 
Knowledge  of  tae  mode.  However,  if  we  use  Rayleigh’s  Principle,  we  can 


INifWjw.wn, ,  *,  >■  ’  "  i  'i'  -  •iti.v.v,  -it'hW 


translate  accurate  Information  about  the  code  into  accurate  information 
about  the  frequency.  This  is  exactly  what  Rayleigh  did  for  spherical  domes 
and  cylinders,  and  w«  shall  derive  general  formulas  for  domes  with  the  edge 
conditions  of  these  two  cases  in  Section 

Case  (III):  n„_  *  N  *w*»  *  0  at  a  *  a0 

ss  s  ss 

This  is  the  "freest”  of  the  remaining  boundary  conditione.  The  analysis 
proceeds  as  in  Case  (x)  except  that  (25)  is  replaced  by 
V  =  ijc(a)  *  B2s(a)  *  B^i*^  ♦  By/  ~  C. 

After  eliainating  the  matrix  equation  of  tha  system  is 

(O2*,,  *  A-‘X„)  1 

2“l/2A2  AX„,  (a2AX^*X2tt  >  B,  =  [0] 


nn  xi  2 

(02X13  ♦  A~kXtk) 

2" 1/2 A2  AX22 

(a2 a  x2j  x2«,  > 

It  2  i-ly 

*  A  Aj2 

(02A'*Xj3  ♦  A’^x,,,) 

where 

n  <*>c<a) 

V  ,  ss 

n_stQ)c(a) 

!!  «„ 

,  •<*> 

•  *2i 

fGn 

n  S<a>s(a) 

v  .  .  ,  ” 

32  • 

fG 

n 

A  33  * 

fGn 

X  S  yf 
2* 


,i(a*) 


and  tha  remaining  X's  are  defined  as  in  Case  (1).  We  find  for  the  frequency 
X  X 

Q*A  *  - -  *  o* 

■  V  * \>u«  •  V 

and  after  some  reduction 
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Q  *  anIA'l/2  *  aiIl{«V«0»lA 

The  frequency  1«  now  higher  by  a  factor  of  roughly  c“1'2  than  in  Case  (II). 
The  eode  is 

w  *  j)  -  xco*  C  »  u  *  u^7^(«)  ,  v  »  v(7>(0) 

bs  *  A(o)^cos(c  ♦  (*/•♦)}  ,  *  b^^'^o)  -  Mxcos  t 


ss 


ee 


-A  2(o)x 


LVJ 


tin  C  ,  "se  *  'A(a)^(i  -v  bluest  C  ♦  (*/**)} 


qs  »  -A3<o)*»in{{  ♦  (*/“)}  *  q  *  A2<<? )M*sir=c 


ss 


a  -A*  (o)Gn(o)x 


ain{(  ♦  (»/**)) 

i>j  /  ltj  r  n  <s>, 

♦  A~1(o0>w'7^(Ofi)Gn(oA 

l 

n  *  -G  xcos  c 

nA 


/ 1 

)J"  i'i  <  L,  j  “ 

i  Ln».  !  (i,i 


♦  L 


A 


where 


X  -  {h(o)/H(o# )  )eS/7^(o0 ) 


o  =  o0 


This  mode  differs  from  those  in  the  two  preceding  cases  in  two  important 
ways.  First,  the  dispiacenents  are  no  longer  completely  inextensional,  for 
the  bending  solutions  make  a  contribution  to  w  neat*  the  edge.  Second  the 
effect  of  the  mesbrane  solution  is  not  now  completely  negligible  but  is 
felt  in  the  formulas  for  the  direct  stresses* 


L«  - 


U  = 


"“ss 


(a)  . 


fn 


(a) 


»e 


n  <Q>n  <*>  -  n  (Q)n  <*> 
SS  S0  SB  SS 


<S> 


^se  5 


Mn 


ss 

(5) 


ss 


<“>n  <s>-.  <>> 

s8 


n  n  „ 
ss  SB 


ss 
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Can  { £¥): 


■S' 

an* 

i 

i. 


& 


& 

9- 

i- 

r 

s< 

y 

f 

Sr 


VC 


Can  <  £V );  »  M#fl  «  w  ■  D  «  0  at  o  »  ot 

Tha  analysis  is  tha  «m  in  this  can  as  in  Can  (III)  with  an  obvious 
chatg*  in  tha  last  boundary  condition.  Tha  natural  fraqucncy  is  found  to 

b* 


Q2  ■*  ajy2!*’1 


(35) 


&iv 


2  *  2^^ 


(7) 


"lil  ft..<s> ' 


<s> 


1 


»  (0’»  .<s>  -  \ 

ss  s8  Sfl  *•  J 


a  -  a. 


and  tha  soda  Is 

•*  «  tt*7*(o)  *xsin(c  -  (s/1*))  »  u  «  u^7*(o)  »  v  *  v*7\0)J 

S>a>  -A(a)y*in  (  •  b#  «  b_*7*(o)  ♦  x  sinfc  -  (*/*»>} 

oos(;  -  (*/*»))  ,  *  A(o)Mx(l-w)sinc 

q^  »  -A*(o)xcoa  t  «  q#  •  A2(o)Mxcos(c  -  (*/•»)} 

xcos  c 

n,s(5)(ol 


r.  *1 

1 

V 

i 

S 

j  71  ~?,  (o))£ 

1 — 

a 

a> 

L®. 

t 

_v. 

C  n _ 

jY 

;  ,  S3 

*  -a"1  VolGy^Cf)  j 

f 

w 

-  (  -t\ . 

(36) 


♦  L„ 


®*8 


<4) 


(a) 


*  jj 


whsra 


n  _  3  G  (o)x»ia(c  -  (*/•*)) 
89  c 


2l/2{H(o)/rf(oft))  a«w(7,io0) 


Vnis  fraquancy  and  soda  ara  quilitativaly  auch  liks  thon  in  Cass  (III). 
We  saa  fro*  (35)  and  (33)  that  tha  fraquancy  estiaat*  in  tha  praswnt  casa  is 
largsr  than  in  Casa  (III)  by  a  siapla  factor  22^2. 


sT 

& 
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The  edge  conditions  considered  in  Cases  (I)  -  (IV)  all  have  ngg  =  Ng(J  =  0 
at  J  1  °0»  ie  the  edges  of  the  shell  have  been  free  to  move  in  directions 
tangent  to  the  middle  surface.  We  have  now  exhausted  all  the  cases  with  this 
property.  If  we  work  out  similar  analyses  for 
Ca-'  (V):  ngs  =  v  =  Qg  =  mgg  =  0  at  a  =  oQ 
Case  (VI):  u  =  N'sB  =  0S  =  mss  -  0  at  o-aQ, 
we  see  tnat  natural  frequencies  in  the  range  J)2  >  0  (X-1)  cannot  occur,  ie.  for 
these  edge  conditions  ail  the  natural  frequencies  obey 

a2  ^  od). 

But  all  tne  remaining  edge  conditions  are  obtained  from  (V)  or  (VI)  by  "tightenin'’" 
some  of  tne  conditions,  ilence  in  all  the  remaining  cases  the  natural  frequencies 
are  at  leas*  as  high  as  in  Cases  (V)  or  (VI).  We  conclude  that  only  for  Cases 
(I)  -  (IV)  can  we  find  natural  frequencies  in  the  range 
n2  1  c(c1/2). 

Now  it  is  easy  to  see  from  the  motion  equations  that  inextensaonal  solutions, 
i,e.  solutions  naving  the  property  that 

nss»  ne3»  nsdL  ^ 

and  all  other  quantities  are  0(1), cannot  occur  when  32  >  0(1).  Hence  ror  a  dome 
inextensional  modes  and  frequencies  can  occur  only  in  Cases  (I)  -  (IV),  only 
wnen  tna  edge  is  free  to  move  tangentially. 
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<4  APPLICATIONS  Of  RAYLEIGH'S  PRINCIPLE 


In  this  Section  we  shall  see  how  estimates  of  the  inextensional  frequencies 
can  be  obtained  for  the  Cases  (I)  and  (II)  by  using  Raleigh's  Principle. 

Rayleigh's  Principle  states  that  for  any  field  of  displacements  satisfying 
the  edge  conditions  on  displacements, 

fi2  1  °£2  S  {Es  ♦  Es)  !  K  (37) 

where  £s,  Eg  and  K  are  to  be  calculated  free,  the  given  field  of  displacements  by 

means  of  (13).  The  accuracy  of  the  estimated  rrequency,  Q_,  depends  on  (and  is 
usually  much  oetter  than)  the  accuracy  of  the  assimed  displacement  field.  -fe 
must  emphasize  (because  it  is  occasionally  overlooked)  tne  effect  of  the  edge 
conditions  on  tne  displacements.  If  these  edge  conditions  are  not  satisfied 
by  the  chosen  displacement  field,  may  differ  wildly  from  Q  and  need  not  even 
be  tne  larger  of  the  two. 

In  applying  Rayleigh's  Principle  to  Cases  (I)  and  (II)  we  shall  take  as 
the  trial  displacements  the  approximate  modes  given  for  these  Cases  by  our  pre¬ 
vious  analysis.  From  (13),  (*♦)  and  (5)  we  have 

Es  =  (i  -  v2)'2/o®o(nss2  ♦  n692  -  2vnssnfte  ♦  2(1  ♦v^e2)  rQsin#do  (38) 
Ea  =  »-(l  -v2)*1  /";0{mss2  ♦  aee2  -  2«nssm8e  ♦  2(1  ♦  vh.^2 )r0sin*dn  (39) 

K  =  /°°  (u2  ♦  v2  ♦  w2)  r.sinfdo  (**-  > 

o=0  0 

Referring  to  the  formulas  (30)  and  (32)  for  tne  modes  in  the  two  cases,  we  see 
that  two  kinds  often*  occur,  namely  terms  of  inextensional  and  edge-effect  types. 
The  integrals  of  the  inextensional  terms  are  of  the  same  order  as  the  tenr.s  then- 
selves  and  cannot  be  evaluated  explicitly  until  the  shell  shape  is  specified. 

The  integrals  of  the  edge  effect  terms  are  smaller  by  an  order  of  magnitude 
thar  the  terms  themselves  and  can  be  evaluated  explicitly  (though  approximately ) 
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by  use  of  the  Laplace  Method  for  asymptotic  approximation  of  definite  integrals. 


for  example  in  Case  (ID,  (33)  and  (32)  lead  to 
Es  *  (.1  -  v2)”2 

a  (1  -  v2)"2  /*•  [A"2  (o)x2(c)  re(a)sin#(o)^,2(o)3  cos2  [c  -  (*/4)]do 
-  /°3{  A"2(o)r  ”‘(o)  sin  *(c  )A  (<j)[D(7)(o0  )]2)  * 
cos2[c  -  (a/**)]  do 


•mere 

;.(i)  =  ii(i  )A(o  )/[h(o0  )A(oc  )J 

T;ie  function  e7 "•  nas  the  value  unity  !nr  o  -  oQand  decreases  ru.  ---/  zero  as 
o  decrease  frcr-,  oQ  .  Hence  tiiis  integral  is  of  Laplace  type,  v*.  jni/  "“e  region 
near  a-  j,  contributes  appreciably  t.  the  integral.  We  "nay  tee : »•  ‘  re  appoxiiate 
it  oy 

L  1  { A~2r  _1n^7,2aintt  v.l/2)/°°  e2^(l  +  sin  2c)  da 

J  0  3 

Tr.e  Integra-  in  this  expression  can  ne  evaluated  approximately  «•'**  •  *.»  aid  or  (31) 

*  ■  i  ve 

T.  >  { A_3r  ‘lO(7)2sinti,  (2,/?/8) 

3  U  Oq 

We  know  t:iat  ,.’2  -  3(x”3),  nence 

Cq  )  1  *r^  1  * '  (oq  )• 

“.us,  ‘ir.aily  wc  find 

L  =  (1  -'8  K'^ir^Og  )]1/2[P*  7  hj0  )]2  sin$(og) 

In  a  similar  manner  may  oe  estimated, 

•  (J/3)a  J[2rg(og)]1/,?[3^  7^(a-j)J2  sin*  (og) 

Hence  for  .  ase  f 11 ) 

r  ’r ra(3g)/2]  i/2[3(  7)(  33 )  j2  sin  *(o0) 

K(?)  (41) 


20 


wnere 


k<7)  -  /C|J  {  [«i<7-(o)]  2  ♦  [v^72(o)]2  •*•  [w*7'(o)]2}  rAsin#do 
o  «0  0 

Applying  the  same  analysis  in  Case  (I)  we  find 

E  =  OU-*) 
s 

EB  *  E3(?)  *  °  <X"5) 

e0(7)  =  rV°*<[»  <7>l2  4  [mJ7>]2  -  (7)=J7) 


b  £So“-ss  J  ‘■“ee  J  ss  89 

♦  2(1  ♦  w)la  a^7^j2)  rfl  sin#  do 


K  =  K' 


Hence  we  obtain  tne  estimate 

a2E  =  Eb(7'Vk(?)  s  0(A'‘f)  (42) 

We  see  that  in  this  case  the  estimate  given  by  Rayleigh's  Principle  can  be 
derived  solely  from  the  irextensional  displacements.  This  is  not  true  of  the 
estimate  just  obtained  in  Case  (II),  nor  is  it  true  of  the  Rayleigh  estimates 
that  are  obtained  for  the  inextensional  frequencies  in  Cases  (III)  and  (IV). 

Equation  (*»2)  is  of  course  just  the  estimate  that  Rayleigh  used  to  'ind  the 
inextensional  frequencies  for  a  spherical  dome.  However,  neither  Rayleigh  nor 
any  subsequent  investigator  seems  to  have  been  sure  of  the  conditions  under  which 
the  estimate  is  accurate.  We  now  see  that  it  is  accurate  only  when  the  edge  of 
the  dome  is  free. 
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S>.  INEXTENSIONAL  f’R£')UENCIES  TOR  A  SfHERICAE  DOffi 

In  this  Section  we  carry  out  the  calculation  of  the  two  lowest  inaxtensionai 
frequencies  for  spnerical  dcKe  under  the  edge  conditions  of  Case*  (II)  and  (III), 
using  formulas  < 1 }  and  (33),  respectively. 


for  a  spherical  dime  the  inextensional  solution  that  is  finite  at  4  *  0  has 


u(7)  -  v(7>  -  s£n  ^tanra  (*/2),  m  ^  2 

w(?)  =  -(n  ▼  cos+)  tan®  (4/ 2) 


(43) 


P^7^  =  }sin$  -  iura  ♦  cos$)  csc$}tan°  ($/2) 


and  tno  sinetic  energy  is  given  by 

K^7*  =  K(*  ,  m)  =  tan2*  ($/2){  2sin2*  ♦  (n  ♦  cosf)2isinf  d#  . 

0  0 

'j/leigr.  has  shown  how  this  integral  can  be  evaluated  for  integer  values  of  m. 

In  Case  (II)  Equation  fal)  reduces  to 

,2  =  2-'/2t3/2ain#0  tan2®  (*0/2)(  sir-40  -x  (m  ♦  cos#())csc  *0}2  (44) 

M  1  —  I  1  ■ 

dK(#c,-i) 

iigare  1  snows  graphs  of  the  relations  between  and  $0  for  m  =  2  and  3,  ob¬ 
tained  fro".  I au). 

In  Case  (III)  the  frequency  is  given  by  (33).  To  evaluate  this  for  a 
spnere,  we  observe  first  that  the  membrane  solution  finite  at  4  =  0  has  (see 

•SOiS  [0  ] ) 


n_  i  ^  *-n  (  , 

s3  ss 


and  (33)  reduces  to 
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we  Must  solve  the  system  of  equations  obtained  from  the  notion  aquations  by 
setting  t  =  0  and  taking  for  u,  v  and  w  the  inextensional  displacements  (43). 
The  governing  equation  is 
dz 

-  +  (m  ♦  2  cos+)  csct  z  1  —(1  -  v2){u  ♦  v  -  w(m  ♦  cost)csct) 

dt 

and  a  particular  solution  (which  is  all  we  need)  is 

z(to»  1  -(1  -  v2)K  (#0,a)  sin-2t3  tan"®  ( t^/2 ) 

Combining  this  with  (45)  and  (33)  we  find 

.T2  3  cl/2(a  cost)2  sin#0  tan2®  f0 
2l/2K(t0,m) 

The  frequencies  predicted  by  (46)  when  ra  =  2  and  3  are  shown  ir.  Tigure  2. 

The  inextensional  frequencies  in  Case  (IV)  are  2,/2tiffi*s  those  of  Case  (II- 
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b.  nisei  ssi ok 


The  results  derived  in  the  preceding  Sections  aid  our  understanding  of 
inextensional  modes  and  the  roles  they  play  in  shell  vibration  problems.  In 
brief,  wo  may  say  that,  when  inextensional  modes  can  occur,  they  are  associated 
with  frequencies  lower  than  those  associated  with  any  other  type  of  mode.  At 
most  one  inextensional  frequency  is  found  for  each  m  >  2,  but  for  a  doio  none 
can  occur  if  the  edge  conditions  involve  significant  constraint  against  no* : on 
tangent  to  tne  shell  surface.  Tne  inextensional  rreqjencies  are  rar  -'ore 
sensitive  to  tne  edge  conditions  than  are  the  (higher)  membrane  and  bending 
frequencies. 

We  have  seen  that  the  procedure  used  by  Rayleigh  to  find  inextensional  rre- 
quencies  yields  very  nearly  tne  results  obtained  by  the  present  metKtd  for  a 
free  edre,  Casa  (I).  Tor  Cases  (II)  -  ( xJ  )  the  present  method  predicts  modes 
with  predominantly  inextensional  displacements,  which  cannot  be  found  hv  'Jayl»i'h*s 
procedure.  Tne  frequencies  in  these  Cases  are  higher  than  in  "ase  (T)  but  still 
i  ;»  oirared  to  the  lowest  frequencies  obtained  for  all  the  rr^ainir'  ed*e 
conditions. 

An  interesting  aspect  of  tnis  analysis  is  this.  When  <  3(i  2),  the 
only  quantities  among  all  tne  eight  solutions  that  depend  om  ,T  are  the  direct 
stresses  associated  witn  the  inextensional  solutions.  When  the-e  smlutions 
are  inserted  in  the  boundary  conditions  and  a  natural  frequency  is  calculated, 
it  is  clear  that  tne  inextensional  direct  stresses  are  the  indispensable  ingred¬ 
ient  ;  of  the  calculation.  Tor,  if  they  are  absent  from  the  frequency  equation, 
it  cme-;  not  contain  tne  frequency  and  cannot  be  satisfied,  and  no  natural 
frequencies  will  be  found.  Yet,  despite  the  importance  of  tnese  direct  stresses, 

2*» 


'  n  ^  '  '*  "j '  *  *  ***-»><  *v 1  ^ ^ 


b 


they  can  be  completely  neglected  in  applying  Rayleigh's  Principle  to  find  the 
inextensional  frequencies  for  a  free-edged  done. 

He  have  not  considered  boundary  conditions  of  elastic  constraint  at  the 
edge.  In  general  we  may  expect  that  these  will  produce  frequencies  lying  betw*«.n 
tnose  associated  with  the  two  "pure"  edge  conditions  that  are  combined  to  give 
tne  elastic  condition,  for  example,  the  lowsst  natural  frequency  associated 
with  the  boundary  condition 

”s»  ■  “S8  ="  =t-s.  *  (1  -«»  ■ 
wnere  0  <  £  <  1,  should  satisfy 

auih~l/2  1  0  I“ivrl/2  • 

Although  we  have  chosen  to  demonstrate  this  procedure  for  don»es,  it  ought 
to  work  equally  well  for  shells  with  two  edges.  However,  it  remains  always 
subject  to  the  condition  thet  »2i  «  1. 


7  5 


-are  i.. 


Inextensicnai  frequencies  as  functions  of  Lape 

&r  tne  t4fe  Conaition  n  =  N  =0=u  =  J 

ss  se  s 
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In  this  paper,  the  non-synetric ,  free,  elastic  vibrations  of  thin  doties  of 
revolution  are  studied.  It  is  assused  that  the  frequency  is  low.  The  asymptotic 
approximations  previously  giveo  by  the  writer  are  used  to  estimate  the  general 
solution  to  the  shell  vibration  equations  at  low  frequencies.  Approximations 
for  the  low  natural  frequencies  and  nodes  are  derived  systematically  under  a 
variety  of  edge  conditions.  Low  natural  frequencies  are  found  only  when  the 
edge  conditions  impose  no  forces  tangent  to  the  shell  surface,  When  the  edge 
is  free  (and  only  then!  Rayleigh’s  inexteosional  frequencies  are  recovered. 

For  certain  other  edge  conditions  new  natural  frequencies  are  found  that  are 
above  Rayleigh's  frequencies  but  still  low  compared,  e.g.,  with  the  lowest 
membrane  frequency.  The  displacement  modes  associated  with  these  new  frequencies 
are  mostly  of  inextensional  type.  The  general  results  are  applied  to  estimate 
these  new  frequencies  for  spherical  douiea. 
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